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The adiabatic theorem provides the basis for the adiabatic model of quantum computation. Re- 
cently the conditions required for the adiabatic theorem to hold have become a subject of some 
controversy. Here we show that the reported violations of the adiabatic theorem all arise from 
resonant transitions between energy levels. In the absence of fast driven oscillations the traditional 
adiabatic theorem holds. Implications for adiabatic quantum computation is discussed. 



A statement of the traditional adiabatic theorem [l|, 0, 
0], as described in most recent pubhcations, is as follows: 
Consider a system with a time dependent Hamiltonian 
H{t) and a wave function \%l}{t)), which is the solution of 
the Schrodinger equation [fi — 1) 



i\m) = H{t)m)). 



(1) 



Let \En{t)) be the instantaneous eigenstates of H{t) with 
eigenvalues En{t). If at an initial time t ~ the sys- 
tem starts in an eigenstate [^^^(O)} of the Hamiltonian 
H{0), it will remain in the same instantaneous eigenstate, 
\En{t)), at a later time t — T, as long as the evolution of 
the Hamiltonian is slow enough to satisfy 



max 

tG[0,T] 



(s,„(t)|£;„(t)} 



-^nm (^) 



< 1 



for all m n, (2) 



where Enm{t) = En{t) — Em{t). One can easily show 
that: {Em{t)\Mt)) = {Emit)\H\E„it))/Er,mit). The 
adiabatic theorem has recently gained renewed attention 
as it provides the basis for one of the important schemes 
for quantum computation, i.e., adiabatic quantum com- 
putation 0, [l]. 

Recently, the adiabatic condition ^ has become a sub- 
ject of controversy. It was first shown by Marzlin and 
Sanders Q and then by Tong et al. 0] that if a first 
system with Hamiltonian H{t) follows an adiabatic evo- 
lution, a second system defined by Hamiltonian 



H{t) = -U\t)H{t)U{t), U{t) EE Te^'/o 



cannot have an adiabatic evolution unless 



|(i?„(0|i?„(0))|«l. 



(3) 



(4) 



even if both systems satisfy the same condition ([2]). Here, 
T denotes time ordering operator. Recently, the validity 
of the adiabatic theorem was experimentally examined 
Q , and ^ was reported to be neither sufficient nor nec- 
essary condition for adiabaticity. 

These inconsistencies have created debates among re- 
searchers [13, 

17l. llSlI , reexamma- 



ternative conditions 
tion of some adiabatic algorithms [19| , or generalizations 
of the adiabatic theorem to open quantum systems [IT 



While it is valuable to find new conditions that guaran- 
tee adiabaticity in general, it is important to understand 
why the traditional adiabatic condition ^ is sufficient 
for some Hamiltonians, but not for others. Moreover, 
from the practical point of view it is much easier to work 
with a simple condition like ([2]) than some other more 
sophisticated ones. In this letter, we relate the reported 
violations of the traditional adiabatic theorem to reso- 
nant transitions between energy levels. We further show 
that in the absence of such resonant oscillations, the tra- 
ditional adiabatic condition is sufficient to assure adia- 
baticity. Our line of thought is close to that of Duki et 
al. 9], but largely extended with rigorous mathematical 
proofs. 

It is well known that fast driven oscillations invalidate 
the adiabatic theorem [^ij. Consider a simple example 
of a two-state system driven by an oscillatory force: 



1 



-eUz — V smujQt 



(5) 



We take V to be a small positive number. The exact 
instantaneous eigenvalues and eigenstates are 



£^0,1 = 



|£^o,i)= 



(6) 



where VL = ^ 41/2 gin^ and a± = ^(r2±e)/217. 
To the lowest order in F, « e -f {IV'^ j uji^t, 
a+ « 1 - (FV2e^) sin^ ujqI, and a" « (T^/e) sin^ot. The 
traditional adiabatic condition Q leads to 



COSWotI < 1, 



(7) 



which is satisfied if Vujq <C e^. Near resonance (o^o ~ e), 
this requires V <C e, t^^o- The adiabatic theorem therefore 
states that if at < = the system starts in its ground 
state |i?o(0)) = (1,0)"^, it will stay in the instantaneous 
ground state at later times. This, however, is not true as 
we shall see below. 

Using the rotating wave approximation, the wave func- 
tion of the system at resonance (e = loq) is given by 



m)) 



-iujQt/2 , 



;yt/2 \ 
?LVt/2 ) 



(8) 
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Therefore, the ground state probabihty 

Po{t) = |(£;o(i)|V^W)l' « (cosyi + l)/2 (9) 

oscillates with the Rabi frequency fn — V/2tt. At time 
T ~ Tji/2 — tt/V, the system will be in the excited state 
with probability Pi = 1, violating the adiabatic theorem. 
Reducing the oscillation amplitude V would only increase 
the Rabi period Tr, and does not keep the system in the 
ground state. Therefore, adiabaticity is only satisfied for 
a time T <^ Tr. Indeed some new versions of adiabatic 
condition set an upper bound on the time T in order to 



guarantee adiabaticity [15|,[16|. However, as we shall see, 
this is not necessary in general. Before that, let us take 
a close look at the inconsistency discussed in [1, 0| • 

Let us assume that H{t) is a slowly varying Hamilto- 
nian for which the adiabatic theorem holds. This means 
that if at time t = 0, the system starts in an eigenstate 
\E°) (= \En{0))) of i?(0), at time t, the wave function of 
the system will be (see below for proof) 



(10) 



Hereafter, we use a gauge in which {En{t)\En{t)) — 0. 

Now consider another system with Hamiltonian ([3]). 
The eigenvalues and eigenstates of the new Hamiltonian 
are E„{t) = -E„{t) and |^„(t)) = W{t)\E„{t)), respec- 
tively. From pn)) . we have 



|^„(i))«e^/o^"(*')'^*'|£;JJ). 



(11) 



It was shown in Refs. 01 that for system H the adi- 
abatic theorem holds only when ^ holds, even if the 
adiabatic condition ([2|) is satisfied. To understand this, 
let us write H in the basis \E^): 



m^Y.^El\H{t)\El)\El){El 



(12) 



does not keep the system in its original eigenstate beyond 
half a Rabi period. Notice that Eq. (|4|) is equivalent to 
\En{t)) K. e"^(*)|i?°), where (t>{t) is some time-dependent 
phase. This leads to {E^{t)\E„{t)) cx {E°JE°) = 0. 
Therefore, the oscillatory terms in (|14p will all vanish if 
dH) is satisfiedjleading to an adiabatic evolution in agree- 
ment with 0] . 

We now provide a general proof for the adiabatic theo- 
rem emphasizing the role of resonant transitions. Let us 
write the wave function of the system as: 



\m) 



^a„(Oe-'/o^"(*')'^*'|£;„(i)). 



(15) 



For a time-independent Hamiltonian, a„(t) is a constant 
while for a slowly varying Hamiltonian it is a slow func- 
tion of time. Substituting p5p into the Schrodinger equa- 
tion ([T]), we find 



dt 



^m{t) = -J2 an{t){E^m{t)\Enit))e'"^o''-^''^ 

Integrating over t, we get 
am{T) - a™(0) = 
-J2 r dt a„(t)(i?™(t)|£;„(t)}e-^/o 

/ -'0 



(16) 



dt 



To assure adiabaticity, the right hand side of this equa- 
tion should be small. With the initial condition am{0) — 
Smn, this would immediately yield (jlOp . Since the expo- 
nential term in the integrand of p6p is a rapidly oscil- 
lating function, if the rest of the terms vary very slowly, 
the integral will be small. To make this statement more 
quantitative, let us define the right hand side of the above 



equation as the error £,„ = — ^ 
batic evolution, where 



£nm for the adia- 



However 



{E"jH{t)\E'J 



Using pO)) we find 



^{E^jU^H{t)U{t)\E'„) 

-i{i^m{t)\ijn{t)). (13) 



H{t) = -J2En{t)\E^,){E'J 



n.m 



{E^it)\MmE"J{E'J. 



(14) 



where ujnm{t) = 7/0 Enm{t')dt' . The second line in 
(|14p involves rapidly oscillating terms that cause reso- 
nant transitions between the levels. The amplitude of 
each oscillatory term is \{Em{t)\En{t))\- Hence satisfy- 
ing ^ will only reduce this amplitude and, as we saw be- 
fore, it does not eliminate Rabi oscillations and therefore 



ciM„™ (*)£;„„ (t)e-^/o^"'"'*')'**', (17) 



and 



{E„MEn{t)) 
Enm (^) 



Using the Fourier transformation: 
dte'^'Anmit), we find 



duj 1 

— J dtAnm{^)Enm{t)e' 



.(t)]t 



(18) 



(19) 



The integral in (|T9)) is suppressed by the oscillatory 
exponential in the integrand, except along a path in 
the two dimensional t-uo plane defined by the equation 
uj = LiJnm{t)i where there is no oscillation. In the presence 
of resonant oscillations, Anra{^) has finite components at 



3 



frequencies LOnmit), hence the contribution from such a 
dominant path becomes 



'o 

< T maJi\Anm{i^nm{t))EnTn{t)\- 



(20) 



The error therefore grows as a function of T. As a resuh, 
to assure adiabaticity one needs an upper Hmit for T, as 
expected for the case of resonant osciUations. However, 
this is not the case for a general system without resonant 
oscillations, as we shall see below. 

In the absence of such oscillations, Anm{t) can be 
made as slow as desired by making the evolution slow. 
In Fourier space, this makes Anmi^) sharply peaked 
at low frequencies with a cutoff frequency ujc propor- 
tional to the rate of change of the Hamiltonian. To see 
this, let us introduce a new variable s — t/T. Since 
an{t) = aTi(O) + 0(£„), to the lowest order in the small 
error £„ we have [22| 



Jo 



. \E^{s)\d/ds\E^{s)) 



(21) 



where Co = luT is the dimensionless frequency. The inte- 
gral on the right hand side is independent of T. Let Coc 
be the largest dimensionless frequency of Anm i^j) ■ There- 
fore LUc = Cbc/T can be made arbitrarily small by making 
T large. Notice that Cjc is a constant that only depends 
on the properties of the Hamiltonian and does not de- 
pend on the evolution time T. 

If ojc LOnm (t) , One can neglect u> in the exponential 
in the integrand of (|19p and therefore perform the t- and 
tj-integrations independently, yielding 



(0)1 <Wc/ \An^{t)\dt 

Jo 



max I A^i^ (^) 



< LOc max 
t 



{E^{t)\En{t)) 



(22) 



Therefore, Snm can be made arbitrarily small by only 
satisfying the adiabatic condition 

The same conclusion can also be reached from a differ- 
ent angle. Using integration by parts, Eq. (fT7| becomes 



^™(r)e-'r^-(*')^*'-A™(0) 



< |^„„(r)|-|-|A„™(0)|-f / dt\Anmit)\. (23) 

Jo 

The last term above is responsible for the breakdown of 
the adiabatic theorem. Let us define i^, i = 1, Mnm, as 



the solutions to Anmiti) = 0, where Mnm is the number 
of zeros of Anm{t) in the interval [0,T]. Since Anm(t) is 
monotonic between ti and ii+i, we can write 



dt\Anrnit)\ = ^ 



i=0 



dtAnm {t) 



\A„m{U+l) - Ar^,Jl{U)\ ,{24) 



i=0 

where we have defined to = and iAf„„+i = T. Thus 

£nT?i 2 > \Anm{ti)\ < 2Mnm max |j4„„j(i)| 
1=0 



< 2 Mnm max 

te[o,T] 



{Em{t)\Ea{t)) 



(25) 



Since the error depends on Mnm, it is important to un- 
derstand how Mnm depends on the evolution time T. 

Let us first consider a Hamiltonian that has an oscilla- 
tory term with frequency loq. Oscillations of the Hamil- 
tonian will create oscillations in Anm{t) and therefore 
the number of zeros of Anm{t) will grow with time as 
Mnm woT. In that case, without putting an upper 
bound on T, it is not possible to limit the error Snm- 
This is in agreement with our previous observation for 
cases involving resonant transitions, as well as the ad- 



T1,[T1]. On the 



ditional conditions introduced in Refs. 
other hand, if Mnm does not grow with time, one can 
always reduce Enm by just satisfying ^ without a need 
to limit T. To see this, let us again use the dimensionless 
parameter s — t/T. If by slowing down the evolution, we 
only change T and not other parameters in the Hamilto- 
nian, then the Hamiltonian H{s) and its eigenvalue En{s) 
and eigenfunctions |£'„(s)) will all be independent of T. 
Again to the lowest order in £„, a„(t) w a„(0) and there- 
fore from (UHl), Anm{t) « an{0){Em{t)\En{t))/Enmit) 

[23 |. The times ti can therefore be obtained by solving 



Anm (t^ — 



a„(0) d {Emis)\d/ds\Eois)) 



T2 ds 



0. (26) 



The number of zeros of this equation, i.e. Mnm, is there- 
fore finite and independent of T. In that case, ((25)) as- 
sures that by just satisfying the adiabatic condition ([2]), 
the error e„„j can be made as small as desired. 

From the above proof it becomes evident that the fol- 
lowing way of stating the adiabatic theorem removes 
all the inconsistences: For a Hamiltonian H{s), where 
s = t/T e [0,1], the evolution of the system starting 
from an eigenstate |i?„(0)) is adiabatic provided that 



T ^ max 

sG[04] 



\{Em\dH/ds\En)\ 

E^ 



for all m^n. (27) 



It should be emphasized that our goal in this paper was 
just to study the sufficiency of the adiabatic condition ^ 



4 



for adiabatic evolution without worrying about the scal- 
ing issue. In other words, we do not discuss dependence 
of the error on the system size. Scaling becomes im- 
portant for determining the performance of an adiabatic 
quantum computer. The exponentially large number of 
states involved in the sum in requires £„,„ to be 
exponentially small in order to keep the sum small. For- 
tunately, this does not put a stringent limitation on the 
time T. To understand this, let us write the error as 



< 



1 



\{Ern\dH/ds\E,,)\ 



(28) 



Typically {Em\dH/ds\En) is exponentially small other- 
wise the curvature of the energy levels 



(fEn 

ds^ 



= 2 



\{E,r,\dH/ds\En) 



Enr. 



{Er 



d^H 



En) 



becomes extremely large due to the sum over exponen- 
tially large number of terms. For the simple example 
of adiabatic Grover search [i^l, it is easy to show that 
{Em\H\Eo) = for m > 1, therefore only first two energy 
levels contribute to the adiabatic evolution. For problems 
with local interactions, the matrix elements typically de- 
cay exponentially with the energy separation between the 
states. This can be checked perturbatively near the be- 
ginning and the end of the evolution (using similar tech- 



niques as in Ref. [2j|). It can also be tested numerically 
for systems with not very large size [25| . Such exponen- 
tial suppression of the matrix elements allows only a few 
energy levels to participate in the calculation of the er- 
ror. Especially, in adiabatic quantum computation when 
the gap between the ground state and the first excited 
state becomes much smaller than other energy separa- 
tions, those two states dominantly determine the error 
of the computation and the evolution time can be deter- 
mined by the minimum gap between those, as has been 
confirmed numerically for up to 20 qubits [l^. More 
work is needed to make these statements mathematically 
rigorous. Moreover, a realistic adiabatic quantum com- 
puter will alwa ys c ou p le to an environment. Therefore, 
other methods 2^, 2y| are necessary to study the evolu- 
tion of such open quantum systems. 

To summarize, we have shown that the inconsistencies 
in the traditional adiabatic theorem reported in the lit- 
erature are all closely related to the fact that for systems 
subject to fast driven oscillations, resonant transitions 
between energy levels cannot be suppressed by just re- 
ducing the amplitude of oscillations, although the adia- 
batic condition ^ can be satisfied. Since the amphtude 
of oscillations determine the Rabi frequency, reducing the 
amplitude would only increase the Rabi period. Within 
a time of the order of half a Rabi period, the system 
will undergo a transition out of its original state. Thus, 
the Rabi period sets an upper limit for the total time 
of the adiabatic evolution. On the other hand, if the 



Hamiltonian of the system does not involve any driven 
oscillations, there is no such mechanism to take the sys- 
tem out of its original state and the traditional adiabatic 
condition is adequate to guarantee adiabaticity. 
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